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We compare the nuclear corrections factors from neutrino deep-inelastic scattering (DIS) 
with the ones coming from a standard analysis of nuclear parton distribution functions 
(nPDF). We focus on a discrepancy between the most precise neutrino DIS data from 
NuTeV and the nuclear PDF coming from the analysis of charged lepton DIS and Drell- 
Yan data. 

1 Introduction 

An indispensable part of any prediction for a process measured at a hadron collider such as 
the LHC are the parton distribution functions (PDFs). Because of the importance of PDFs, 
many groups perform and update global analyses of PDFs for protons [U El [3] and for nuclei 
[H O [S] • Proton PDF are determined from data taken not only on protons but from some data 
taken on nuclear targets, mainly deuterium but also heavy nuclei such as lead and iron in case 
of neutrino DIS. Neutrino DIS data is sensitive to the strange quark content of the proton and 
complements newly available LHC data from W- or Z-boson production. 

In order to include the neutrino DIS data in a global fit to help constrain the proton PDF, we 
have to apply a nuclear correction factor. The nuclear correction factor can be obtained either 
from a specific model of nuclear interactions or from an analysis of nuclear parton distribution 
functions (NPDF) based on experimental data. 

Here, we discuss a compatibility of neutrino DIS data with the nuclear correction factors 
obtained from NPDF analysis focusing on the neutrino DIS data from the NuTeV experiment. 

2 Nuclear correction factors from nuclear PDF 

Nuclear correction factors are in general defined as a ratio of an observable in a nuclear process 
and the same observable in a process involving protons. In the following, we discuss two nuclear 
correction factors both related either to the F 2 structure function in neutrino DIS 



R» cc (F 2 ;x,Q 2 ) 



d A +u A + ... 



(1) 



d A,0 + U A,0 + _ _ _ 



or to the F 2 structure function in charged lepton DIS 




(2) 
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The superscript '0' stands for using the free nucleon PDFs ff' n (x, Q) as given below in Eq. (J5J) . 

Nuclear correction factors such as those defined by Eqs. |l|2| can be either extracted from the 
data or calculated using the extracted parton distribution functions. Here we use the nuclear 
PDF from [7] and [5] where the parameterizations of the nuclear parton distributions of partons 
in bound protons at the input scale of Qo = 1.3 GeV are 

xf k {x,Q ) = c Q x^{l-x) c *e c * x {l + e c *x) c ^ (3) 
d(x,Q )/u(x,Q a ) = c x Cl (l~x) C2 + (1 + c 3 x)(l - x) c * , 

where fk — u v , d v , g,u + d, s, s and u, d are a generalization of the parton parameterizations in 
free protons used in the CTEQ proton analysis [Hj. To account for different nuclear targets, the 
coefficients c k are made to be functions of the nucleon number A 

c k -+ c k (A) = c kfi + c fc ,i (1 - A- Ck -*) , k = {1, . . . , 5} . (4) 

From the input distributions, we can construct the PDFs for a general (A, Z)-nucleus 

#' Z) (*,Q) = j f! /A M) + f: /A (*,Q), (5) 

where we relate the distributions of a bound neutron, f^ A (x, Q), to those of a proton by isospin 
symmetry. 

In the analysis, the same standard kinematic cuts Q > 2 GeV and W > 3.5 GeV were applied 
as in jS] and we obtain a fit with x 2 /dof of 0.946 to 708 data points with 32 free parameters 
(for further details see [7])- In Fig. [I] (solid line), we show how the result of our global analysis 




Figure 1: Nuclear correction factors i?^(F2j x, Q 2 ) and R^ C (F2; x, Q 2 ) from the global NPDF 
analysis compared with the corresponding data for iron target. 

of charged lepton data translates into nuclear correction factors and how the nuclear correction 
factors compare with experimental data. As first observed in [TD], the R£, C {F2 \ x, Q 2 ) correction 
factor calculated using Eq. [T] with parton densities from the fit to the charged lepton nuclear 
data, does not describe the NuTeV data well which raises the question if including neutrino DIS 
data in the global analysis corrects this behavior without spoiling the R < ^j(F2; x,Q 2 ) correction 
factor which fits the charged lepton DIS and DY data well. 
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w 


xf±A (M) 


xIa (M) 


total x 2 (M) 





638 (0.90) 




638 (0.90) 


1/7 


645 (0.91) 


4710 (1.50) 


5355 (1.39) 


1/2 


680 (0.96) 


4405 (1.40) 


5085 (1.32) 


1 


736 (1.04) 


4277 (1.36) 


5014 (1.30) 


oo 




4192 (1.33) 


4192 (1.33) 



Table 1: Summary table of a family of compromise fits. 



3 Neutrino DIS 

To analyze the possible discrepancy between the nuclear correction factor R v cc (Fi \ x,Q 2 ) from 
the fit to charged lepton data and the neutrino DIS data, we have included the NuTeV and 
Chorus neutrino DIS cross-section data in the global fit. The 3134 neutrino DIS cross-section 
data points would clearly dominate 708 charged lepton data in the global fit. That is why, we 
have introduced the weight to the neutrino data and set up a series of fits in order to find a 
compromise fit. % 2 /dof for each compromise fit with a different weight of neutrino DIS data is 
listed in Tab. [I] Each global fit with a different weight results in a different nuclear correction 
factor and in Fig. [2] we see that the weight is a suitable parameter which interpolates between 
the fit using only charged lepton data and the fit using only neutrino data (for further details 
see [5]). In order to decide on how well the compromise fits describe the data we use the x 2 




Figure 2: Nuclear correction factors R^/^(F2;x,Q 2 ) and Rq C (F2;x,Q 2 ) for all fits in Tab.[l] 

goodness-of-fit criterion used in [TTJ [5] . We consider a fit a good compromise if its x 2 f° r both 
data subsets, the charged lepton DIS and DY data and the neutrino DIS data, is within 90% 
confidence level of the fits to only charged lepton or neutrino data. 

We define the 90% percentile £90 used to define the 90% confidence level, by 

/•Coo 

/ P( X 2 ,N)d X 2 = 0.90, (6) 
Jo 

where N is the number of degrees of freedom and P(x > N) = 2"/ 2 r(jv/2) — * s ^ ne probability 
distribution. We can assign a 90% confidence level error band to the x 2 of the fits to the charged 
lepton DIS and DY data and to the neutrino DIS data. By looking at the overall x 2 /dof values 
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and at the plots in Fig. [2] one might conclude that the global fit with w = 1/2 describes both 
data well to constitute a compromise fit. The conclusion changes however when we inspect 
separate contributions to the x 2 /dof from different experiments. The change in global \ 2 is 
mostly due to change in x 2 of the DIS scattering on iron which makes all the compromise fits 
incompatible. The conclusion about incompatibility of the fits rests on NuTeV data having 





Figure 3: Nuclear correction factors R e j$~,{F2]X, Q ) and Rq C (F2;x,Q 2 ) where neutrino data 
were included with uncorrelated systematic errors. 

small errors which can be demonstrated by neglecting the correlations in systematic errors 
which results in a compatible fit of all the data (see Fig. [3]). 



4 Conclusion 

A thorough global NPDF analysis of the combined charged lepton and neutrino data leads us 
to conclude that there is no good compromise description of both the data sets simultaneously. 
The differences can be seen in the low and intermediate x regions where the neutrino DIS 
(NuTeV) do not show a strong shadowing effect as the charged lepton data do. The inability 
to describe all data by one consistent framework poses problems for including the NuTeV data 
to the proton PDF analysis. 
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